Astronomy & Astrophysics manuscript no. master 


February 2, 2008 


(DOI: will be inserted by hand later) 





Altair's inclination from line profile analysis'^ 

A. Reiners^'^ and F. Royer^'^ 

^ Astronomy Department, University of California, Berkeley, CA 94720 

e-mail: areinersOastron . berkeley . edu 
^ Hamburger Sternwarte, Universitat Hamburg, Gojenbergsweg 112, 21029 Hamburg, Germany 
^ Observatoire de Geneve, 51 chemin des Maillettes, 1290 Sauverny, Switzerland 

e-mail: f rederic . royerOobs . unige . ch 
* GEPI, CNRS UMR 8111, Observatoire de Paris, 5 place Janssen, 92195 Meudon cedex, France 

February 2, 2008 

Abstract. We present a detailed spectroscopic study of line broadening in the A7IV-V star Aitair. In a wavelength 
region covering 690A we reconstruct the overall broadening profile taking into account more than 650 spectral 
lines. From the broadening profile we determine the projected rotational velocity v sini, derive an upper limit 
for the equatorial velocity v from the shape of the profile and search for signatures of differential rotation. 
Our redetermined value of v sinz is (227 ± 11) kms~^. Measuring the first two zeros of the Fourier transformed 
broadening profile yield no signatures of differential rotation. We derive that Aitair is seen under an inclination 
angle higher than i — 68° and it rotates at i; < 245kms~^ or slower than 53% of breakup velocity on a la level. 
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1. Introduction in such objects effective temperature varies with grav- 

ity over the stell ar surface: this is the gravity darkening 
As the 12th brightest star in the sky, Aitair (a Aql, jyon Zeipellll924|) 

HR 7557, HD 187642) has been extensive ly studied in c j.- • • iu i, r • >• • > i r 

' /I— : -1 I 1'', „ bcrutimzm g the shape ot individual lines, 

the literature. Recently IBuzasi et al.l 11200411 list most of 17717 7 — Tl fTTTl j , • j .i, • r i- 

.' — ; — : — : ' ' — 7^ . . IGulliver et alJ H1994|) determined the inclination an- 

its physical parameters in their introduction. Aitair is an , r r , -j^t i 

. 1 r S^^ V ega from spectroscopic data. In a more general 

A7IV-V main sequence star known to be a fast rotator , , |7~: i ,onno\ u j-u a- -t t • ^■ 4.- 

, .... . r, 1 -1 context, I RemersI L2003) shows the etrects ot inclination 
with a proiected rotational velocity wsmi — 217 kms , i i • t m r r ^ ^ ^ 

rr-^ , . , . . . and gravity darkening on line prohles oi fast rotators 

llRover et al. 20021. At a distance of 5.14 pc lESA. 1997.1. it • t^, • -, u t, j ^i, ^ • ^ r 

' ' ~ ' u — „ 1, 1 I m t'ls I'ourier domain. He showed that m spectra oi 

is one or the nearest bright stars, allowmg lvan Belle et al. I r . j. 4. 4.1, • i j i i- 1 1 -i 

I 1 ° ' °' ' last rotators the unprojected rotational velocity v is 

(12001') to measure the oblateness of the star using inter- ui j- ^ r u Ti ■ j-u 

' ' . . 11-1 • 1 • /, measurable directly irom the line shape. It is the aim 

ferometric observations. They derived an axial ratio a/b — r ■ , , . 

„„„ , , ,. , . 1 01 this paper to put stronger constraints m deriving the 

1. 140 ± 0.029 and a value for the w sin i of 210 ±13 km s"^ . / r • r , , , 

inclination ot Aitair from spectroscopic measurements by 
Interferometry also allows to make constraints on stellar . , ■ ^ ^\ ^ ■ ^ ^ 1 

. . .1 , , r measuring its (unprojected) equatorial velocity and the 

inclination - a parameter usually not available Irom stan- . . 

. . . . . V sva.1. 

dard spectroscopic techniques. Their determination of the 

inclination i of the rotational axis ranges from 30° to 90° 

within la ( van Belle et al. 20plj Fig. 5). 2 Data 

Spectroscopic techniques are usually not sensitive for 

inclination in "classical" case of spherical stars. On the 2.1. Aitair data 

contrary, for fast rotators, the stellar surface is distorted ai, - i 1 1 j 1. 1 1. nnno 

, .V , ,. , m 1 11 1 Aitair has been observed on Octob er 1st 200 3 with the 

by centrifugal forces and hue profiles depend on the angle j^^ODIE echelle spectrograph iBaranne et al.. .199dl at 

between the rotational axis and the line 01 sight. Moreover, „, ^ . 1 tt j_ ' , 

Observatoire de Haute-Provence. Ihe wavelength cover- 

Send offprint requests to: A. Reiners ^ge spans from 3850 to 6800 A, with a resolving power 

* Based on observations carried out at Observatoire de ^2 000. Five single exposures have been taken, details 

Haute-Provence (OHP), at the German-Spanish Astronomical about the observational data are listed in Table [T] 

Centre, Calar Affo, and retrieved from the Elodie Archive at The orders of the ec helle spectrum have been m erged 

Observatoire de Haute-Provence (OHP) using the method by lErspamer fc Northl l)2002() . The 
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barycentric radial velocity correction varies of about 
0.025 kms^^ in the short time interval of the observa- 
tions, and meanwhile, the peculiar radial velocity of Altair 
is not expected to change significantly. The comparison of 
wavelength positions for telluric lines (O2) insures that 
no significant instrumental wavelength shift occurred be- 
tween the consecutive exposures, and the merged spectra 
are therefore simply co-added. 

2.2. Calibration data 

We used spectra of three stars of similar spectral type to 
calibrate our measurements in Altair's spectrum. Spectra 
observed with the FOCES spectrograph at the 2.2m tele- 
scope at the German-Spanish Astronomical Centre, Calar 
Alto, were taken in February 2002, fo r details on the re- 
duction cf. Reiners fc SchmittI l)2003|) . The spectrum of 
HD 27819 was taken from the ELODIE Archived The 
calibration stars are listed in Tabled 



3. Continuum normalization 

Analyzing spectral line shapes severely depends on the 
determination of the continuum flux. In spectra of stars 
with rotation rates as high as Altair's, spectral lines are 
severely blended and the flux hardly reaches the contin- 
uum in any spectral region. Thus, normalizing Altair's 
spectrum by eye-guided continuum fitting would always 
assume too low a continuum which significantly affects 
the shape of the line profiles. 

Thus we took high quality spectra of other stars of 
similar spectral types, some of them having considerably 
lower rotation rates; in their spectra the continuum is 
nicely visible. From the FOCES campaign we took two 
stars; HD 118623 is an A7III with a rotational velocity 
very similar to Altair. 7 Boo has a rotational velocity of 
V sin I ^ 130kms~^, i.e. significantly lower than that of 
Altair and HD 118 623. In the spectrum of 7 Boo flux 
reaches continuum level in some regions and continuum 
fitting is possible. 

The spectrum of 7 Boo was normalized by fitting the 
regions where flux reaches continuum level. We also used 
this spline fit to normalize the spectrum of HD 118623 
which was taken during the same observing campaign in 
the same instrument setup. The coincidence of their un- 
normalized spectra is shown in Fig.^ giving reason for 
using identical normalization functions. 

With a normalized spectrum of HD 118 623, a star be- 
ing similar to Altair with respect to spectral type and 
projected rotational velocity, we could then normalize the 
spectrum of Altair by matching its overall spectral shape 
to that of HD 118 623. Both normalized spectra are shown 
in Fig[21in the whole wavelength region we used in the fol- 
lowing. 



^ http : //atlas . obs-hp. fr/elodie/ 
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Fig. 1. Small region of the unnormalized spectra of 7 Boo 
(grey) and HD 118 623 (black). The fit to the continuum 
assumed is overplotted and is in good agreement with both 
spectra. 
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Fig. 2. Normalized spectra of HD 118 623 (grey) and 
Altair (black). The whole regions used to derive Altair's 
broadening function is shown. 

4. Reconstructing the global broadening function 

Spectra of stars rotating as fast as Altair are completely 
dominated by rotational broadening. Temperature broad- 
ening and turbulent motions are about two orders of mag- 
nitude smaller than rotation. Rotation affects all spectral 
lines in the same characteristic manner. However, the as- 
sumption we have to make is that the line profiles of in- 
dividual lines do not significantly (at least not systemat- 
ically) vary across the stellar disk. For individual lines 
this is not true, but using hundreds of lines of differ- 
ent species, variations among different lines are expected 
to level out. What one gets is a profile which refiects 
the mean rotational broadening affecting all lines. This 
method is in contrast to studying individual lines (e.g., 
iGulliver et al Il994) and line-specific variations with tem- 
perature or gravity are not taken into account. For the 
case of Altair, systematic effects are not expected to be 
as severe as, e.g., for the case of Vega. We will discuss 
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Table 1. The collected spectra are listed for each stellar target, giving the corresponding spectrograph, the resolving 
power, the Julian date (HJD) and signal-to-no ise ratio (S/N) for each observation. The S/N for ELODIE spectra is 
estimated at 5550 A. Values of v sini are from lRover et all l)2002(l given in kms~^. 



Object 


Sp.Type 


Instrument 


A/AA 


HJD 


S/N 


V sin i 


Altair 


A7IV/V 


ELODIE 


42 000 


2452914.2574 
2452914.2600 
2452914.2628 
2452914.2655 
2452914.2683 


180 
200 
280 
230 
250 


~ 220 


HD 27 819 


A7IV/V 


ELODIE 


42 000 


2450499.4221 


180 


~ 40 


7 Boo 


A7III/IV 


FOCES 


40 000 


2452336.0736 


300 


~ 130 


HD 118 623 


A7III 


FOCES 


40 000 


2452332.1416 


300 


~ 210 



that point in Sect. 15. ll However, individual line profiles 
are severely blended in Altair's spectrum making studies 
of individual lines difficult and wc think that recovering 
line broadening from as many lines as possible yields the 
most reliable information about stellar rotation. 

The reconstruction of Altair's mean broadening pro- 
file is done by a deconvolution algorithm. The algorithm 
searches the function that gives the best fit to the data 
when being convolved with a given template. Thus the 
template is required to be a precise spectrum of the star 
without rotation (neglecting other broadeners). Since no 
such template is available we usually try to simultaneously 
fit the line strengths, i.e., equivalent widths, of all spectral 
lines (cp. Reiners & Schmitt 2003). In the case of Altair, 
however, blending of the large amount of spectral lines 
due to high rotational velocity is so severe, that simulta- 
neously fitting the broadening function and line strengths 
is not possible. We thus did not determine the template 
from the spectrum of Altair but from a spectrum of a star 
of similar spectral type. 

In spectra of slower rotators with less severe blend- 
ing, simultaneously fitting the broadening function and 
line strengths is achievable. From the ELODIE- Archive we 
took a spectrum of HD 27819, a A7IV-V star rotating at 
V sini ^ 40kms~^. Normalization was made by eye using 
the clearly visible continuum. A comparison of the normal- 
ized spectra of Altair and HD 27819 is shown in Fig.|21 
spectral similarity is clearly visible. For this "slow" rota- 
tor HD 27 819 the broadening function and line strengths 
were derived iteratively fro m the spectrum using ou r de- 
convolution algorithm fcp. iReiners fc SchmittlEoO.'^ . We 
used the whole wavelength region between A = 4960 A and 
A — 5650 A. In this spectral region all lines are dominated 
by rotation, i.e., hydrogen lines are not contained. Roughly 
650 spectral lines were used in the fit. Wavelength infor- 
mation was taken f rom the Vienna Atomic Line Database 
(iKunka et al.lll999i VALD). 

The result of this procedure is shown in Fig. 0] where 
we plot the result of convolving the recovered broadening 
function - shown in the inset - with the optimized tem- 
plate. Note that each pixel of the broadening function is a 
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Fig. 3. Comparison of the normalized spectra of 
HD 27819 (grey line) and Altair (black line) in a small 
spectral window. 



free parameter in the algorithm and that the broadening 
function resembles nicely symmetric rotational broaden- 
ing (from the profile we derive a v sini of 46 km s~^). The 
data (grey line) is reproduced in great detail; since indi- 
vidual lines are not severely blended ambiguities between 
the shape of the broadening function and line strengths 
are very small. 

For the deconvolution of Altair's broadening profile we 
now used the same optimized template, i.e., the same line 
strengths as derived for HD 27819. Note that the observed 
spectrum of HD 27819 could not be used as a template 
since it already is broadened by 46kms~^ and a convolu- 
tion of two rotational broadening functions does not yield 
a third one. From the normalized spectrum of Altair we fi- 
nally derived the profile that gives the smallest chi-square 
when convolved with this delta-like template. 

A small region of the result is shown in Fig.|31 fit qual- 
ity is representative for the whole wavelength region. The 
derived broadening function is shown in the inset. Note 
that no additional changes were made to the continuum 
and that the broadening function is nicely symmetric. A 
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Fig. 4. Data (grey) and fit (black) achieved from our dcconvolution algorithm for the case of HD 27819. Fit quality 
is typical for the whole wavelength region. The deconvolved "global" profile is plotted in the inset. 




Fig. 5. Data (grey line) and achieved fit (black line) of Altair. Fit quality is typical for the whole wavelength region. 
The deconvolved "overall" profile is plotted in the inset. 



small linear slope in the overall appearance might be visi- 
ble. While this slope is likely due do tiny continuum mis- 
match, however, it does not significantly affect the line 
shape. The pixel size used is only a fraction of the instru- 
ment's resolution and the small jitter in the broadening 
function is numerical. Taking the jitter as signal-to-noise 
ratio (S/N) yields values of several thousands. We calcu- 
lated the variance of each pixel in the chi-square fit with 
regard to the two nearest neighbors of each pixel. This 
also yields a S/N of several thousands. Unfortunately, un- 
certainties in the broadening function are dominated by 
our knowledge of the continuum and by systematically in- 
correct line strengths (in detail, line strengths of Altair 
definitely differ to those of HD 27819). For the calcula- 
tion in what follows we estimate a rather conservative S/N 
of 1000 as the "Icr" -level. Consequently we use a S/N of 
500 for our calculations of 2a borders. In the left panel of 



Fig. El the global broadening function of Altair is plotted 
with error bars according to a S/N of 1000. 



The global broadening function shows small asymme- 
tries. From our analysis we can not say whether these 
asymmetries are real or due to incorrect line strengths in 
the template. However, within the error bars the broad- 
ening function is symmetric, i.e. the broadening function 
mirrored at the center falls withing the error bars. 

In the right panel of Fig.Elthe Fourier transform of the 
derived broadening function is shown with the respective 
Icr error bars. The signal can clearly be followed up to the 
third sidelobe before it reaches noise level. In the calcu- 
lation of V sin I and i only the positions of the first two 
zeros are needed; apparently both are well defined. 
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Fig. 6. Left panel: global broadening function of Altair reconstructed from about 650 lines. Error bars according to a 
signal-to-noise ratio of 1000 are overplotted. Right panel: Fourier transform of the reconstructed broadening function 
(left panel) with its error bars. Note the logarithmic scale on the y-axis. 



5. Spectroscopic rotation and inclination 

Studying the characteristics of rotational broadening, ze- 
ros of the broadening pr ofile's Fourie r transform are con- 
venient observables (cp. iGravl Il99^ . Since convolutions 
in Fourier domain become multiplications, especially ze- 
ros inferred from any broadening mechanism remain un- 
changed by other effects like turbulence etc.; the intensity 
is changed - the position of the zeros is not. In case of 
a star rotating as fast as Altair, no effects are expected 
inferring signatures at frequencies as low as those inferred 
by rotation. Thus, measuring Fourier transform zeros is a 
very robust method to access rotational broadening - and 
the zeros can easily and precisely be determined. The first 
zero is frequently used to determine stellar projec ted rota- 
tional velocities v sini, as shown in lDravins et al. 1,1990 ) 
for the case of linear limb darkening. 

Effects of rotation and inclination on absorption line 
profiles have been studied by Reiners, (2003). He shows 
that the ratio (72 /^i of the first two zeros qi and (72 is 
sensitive for very fast rotation. He concludes that the ratio 
Q2/Q1 does not depend on the value of v sini, but on the 
equatorial velocity v regardless of the inclination angle i 
(and on the rotation law, see Sect. l5.2|l . Thus measuring 
V sin i from qi and v from the ratio (72 / qi yields the value 
of the inclination angle i\ 

By modeling absorption profiles of fast rotators, 
iReiners (2003.) finds that the ratio 92/91 can be expressed 
as 

92/91 = 1-75 + + (1) 

Values of a and b are given for different stellar models 
in Table 1 of that paper. The mechanism determining the 
values of a and b is gravity darkening described by the 
parameter (3 in 

TeS (X (2) 

with g the stellar surface gravity. For Altair we get a value 
of /3 ~ 0.09 from the tables given in iClaretl H1998i) . i.e. 



the results calculated bv lReiners fc Schmitd (|2003) for the 
case of an FO star are directly applicable to Altair. Thus, 
for the parameters in Eq.^wc find a = .740 • 10^"* and 
b = -.345- 10-6. 

Our measurements of the first two zeros 91 and 92 
of the broadening profile's Fourier transform are given 
in Table|21 Assuming a linear limb darkening law with a 
limb darkening parameter e ~ 0.6 we calculate Altair's 
projected rotational velocity as v sini = 227 ± 2kms-^. 
Note that the formal error due to the uncertainty of 91 
is given. We estimate the systematic error in v sini to be 
about 5 % due to continuum and line strength uncertain- 
ties fcD. iReiners fc SchmittlEool . From Eq.n]we can cal- 
culate Altair's equatorial velocity v from the ratio 92/91- 

In the left panel of Fig.[71 v is plotted as a function 
of 92 /91 with our measurement of v sin i overplotted as a 
horizontal line, and the measurement of 92/91 overplotted 
with its la and 2a error levels as vertical lines. The plot 
shows, that 92/91 is generally expected being smaller than 
1.765, i.e. the measured value of 92/91 = 1-77 is not ex- 
pected for any possible inclination angle. However, within 
the la region we find possible values of ^(92/91) which are 
also consistent with the measured value of v sini. 

Calculating Altair's inclination angle i from the mea- 
sured V sini and w (92/91) we find 1(92/91). We plot 
1(92/91) in the right panel of Fig.[7| la and 2a confidence 
levels in 92/91 and systematic errors of 5% in v sini are 
indicated by dashed lines. We find that our measurements 
of V sini and 92/91 are consistent within the la level for 
values of i > 68°. 

5.1. Intrinsic line variations 

Variations of line profiles and line strengths have not been 
taken into account in our procedure of l ine profile re- 
constr uction. For the case of Vega rAOVL lCxulliver et alJ 
showed that for the two lines they used line intrinsic 
effects dominate the line profile; they observe a flattened 
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Fig. 7. Left panel: Equatorial velocity w as a function of g2/<Zi according to Eq.^ Vertical lines indicate measured 
values of (72/91 and its confidence levels, the horizontal lines shows the measured v sini. The allowed region where 
V > V sini is plotted as a thick line. Right panel: Inclination angle i (in degree) of Altair as function of (72/91 (solid 
thick line). Measurements of (72/(71 are plotted as in left panel. Dashed lines indicate the region allowed taking into 
account the systematic errors of v sini. 



profile instead of a cuspy shape. The difference to our 
method is twofold: (i) gravity darkening as expressed in 
terms of the parameter (3 (Eq.[21) is only 0.09 for our case 
of Altair while it is 0.25 for the case of Vega. Temperature 
variations with gravity are thus much weaker in the case 
of Altair; (ii) our reconstruction procedure uses about 650 
lines to determine the shape of the line profile; we get a 
weighted mean of 650 line profiles. Although we do not 
know specific temperature dependencies for all lines, we 
do not expect them to be systematic. Studying tempera- 
ture and gravity dependence of line strengths and shapes 
for different spectral types would significantly reduce the 
systematic errors but is beyond the scope of this paper. 
For this study we estimate the systematic errors as given 
below and leave the detailed study for a forthcoming pub- 
lication. 



5.2. Differential rotation? 

It has been argued that Altair might undergo substantial 
surface differe ntial rotation since it was found to be an 
X-ray source (|Ferrero et al.lll995h . The parameter (72/(71 
is also a sensitive indicator for solar-like surface differen- 
tial rotation with the equator rotating at higher velocity 
than the polar regions do. In the case of differential ro- 
tation, (72/(71 is expected to be significantly smaller that 
1.75 which is not the case for Altair. We thus conclude 
that no signatures of solar-like differential rotation could 
be found in the line profiles of Altair, i.e. solar-like differ - 
ential is not stronger than 10% ijReiners fc Schm itt"2002^ . 
However, Altair is at least two orders of magnitude fainter 
in X-rays than active F-stars are. Even if a magnetic dy- 
namo was the source for this faint X-ray flux, differential 
rotation may be below the detection limit. 



Table 2. Measured values of Fourier transform zeros and 
calculated parameters v sini, v and i. For v sini formal 
and estimated systematic errors are given. For la and 2a 
confidence levels in v and i, signal-to- noise ratios of the 
broadening function of 1000 and 500 are adopted, respec- 
tively. 



Parameter 


Value 


91 


(2.91 ±0.03) ■ 10"^skm-i 


92 


(5.16 ±0.08) • lO'^skm"^ 


32. 


1.77 ±0.03 


11 




V sini 


(227±2±ll)kms"^ 


V 


< 245kms"^ (la) 




< 305kms-^ (2cr) 


i 


> 68° (la) 




> 45° (2cr) 



6. Conclusions 

We have spectroscopically determined constraints on 
Altair's inclination angle and differential rotation from the 
global rotational broadening profile derived from about 
650 spectral absorption lines. The measured observables 
and our results are summarized in Table|21 where values 
of equatorial velocity v and inclination angle i are given 
with their respective la and 2a levels. 

No signatures of surface differential rotation were de- 
tected, i.e. solar-like diffe rential rotation can be ru led out 
to be stronger than 10% l|Reiners fc SchmittlEoO^ . 

From the measured zeros of the Fourier transformed 
broadening profile we conclude that Altair's equatorial 
velocity is not faster than 245 km s"^ (305 km s^^) or 
53% (72%) of breakup velocity Wprit = 426 ± 12kms-i 
l)van Belle et al. ■■200L) . Thus its rotation axis is seen un- 
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der an inclination angle larger than 68° on a Icr-level (45° 
on a 2a- level). These results put additional constraints on 
the possible rotation rate and inclinatio n angle of Altair 
as de termined from interferometric data ( van Belle et al. | 
l20nil) . 
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